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We apply the recently developed Log N-Log S test of compact star cooling theories for the first 
time to hybrid stars with a color superconducting quark matter core. While there is not yet a 
microscopically founded superconducting quark matter phase which would fulfill constraints from 
cooling phenomenology, we explore the hypothetical 2SC+X phase and show that the magnitude 
and density-dependence of the X-gap can be chosen to satisfy a set of tests: temperature - age (T-t), 
the brightness constraint, Log N-Log S, and the mass spectrum constraint. The latter test appears 
as a new conjecture from the present investigation. 

PACS numbers: 12.38.Mh, 26.60.+C, 95.85.Nv, 97.60. Jd 


I. INTRODUCTION 

Recently, preparations for terrestrial laboratory exper¬ 
iments with heavy-ion collisons have been started, where 
it is planned to access the high-density/ low temperature 
region of the QCD phase diagram and explore physics at 
the phase boundary between hadronic and quark matter, 
e.g., within the CBM experiment at FAIR Darmstadt. 

Predictions for critical parameters in this domain of the 
temperature-density plane are uncertain since they can¬ 
not be checked against Lattice-QCD simulations which 
became rather precise at zero baryon densities. Chi¬ 
ral quark models have been developed and calibrated 
with these results. They can be extended into the finite- 
density domain and suggest a rich structure of color su¬ 
perconducting phases. These hypothetical phase struc¬ 
tures shall imply consequences for the structure and evo¬ 
lution of compact stars, where the constraints from mass 
and radius measurements as well as from the cooling phe¬ 
nomenology have recently reached an unprecedented level 
of precision which allows to develop decisive tests of mod¬ 
els for high-density QCD matter. 

Among compact stars (we will address them also with 
the general term neutron stars (NSs)) one can distigu- 
ish three main classes according to their composition: 
hadron stars, quark stars (bare surfaces or with thin 
crusts), and hybrid stars (HyS). The latter are the sub¬ 
ject of the present study. 

Observations of the surface thermal emission of NSs is 
one of the most promising ways to derive detailed infor¬ 
mation about processes in interiors of compact objects 
(see dli for recent reviews). In |d] (Paper I here¬ 
after) we proposed to use a population synthesis of close- 
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by cooling NSs as an addtional test for theoretical cooling 
curves. This tool, based on calculation of the Log N-Log S 
distribution, was shown to be an effective supplement to 
the standard T-t (Temperature vs. age) test. In Paper 
I we used cooling curves for hadron stars calculated in 
(2]. Here we study cooling curves of HyS calculated in Q 
(Paper II hereafter). 

Except T-t and Log N-Log S we use also the bright¬ 
ness constraint test (BC) suggested in Q. We apply 
altogether three tests - T-t, Log N-Log S, and BC - to 
five sets of cooling curves of HyS. In the next section we 
describe calculation of these curves. In Section III we 
discuss the population synthesis scenario. After that we 
present our results which imply the conjecture of a new 
mass spectrum constraint from Vela-like objects. In Sec¬ 
tion 5 we discuss the results and present our conclusions 
in Section 6. 


II. COOLING CURVES FOR HYBRID STARS 
A. Hybrid stars 

The description of compact star cooling with color su¬ 
perconducting quark matter interior is based on the ap¬ 
proach introduced in Paper II which will be briefly re¬ 
viewed here, see also for a recent summary. A non¬ 
local, chiral quark model is employed which supports 
compact star configurations with a rather large quark 
core due to the relatively low critical densities for the 
deconfinement phase transition from hadronic matter to 
color superconducting quark matter. In the interior of 
the compact star in late cooling stages, when the tem¬ 
perature is well below the opacity temperature T apac ~ 1 
MeV for neutrino untrapping, four phases of quark mat¬ 
ter are possible: normal quark matter (NQ), two-flavor 
superconducting matter (2SC), a mixed phase of both 
(NQ-2SC) and the color-flavor-locking phase (CFL). The 
state-of-the-art calculations for a three-flavor quark mat- 
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ter phase diagram within a chiral (NJL) quark model 
of quark matter and selfconsistently determined quark 
masses are described in Refs. Fill fill Fl3l|. 

The detailed structure of the phase diagram in these 
models still depends on the strength parameter Gd of 
the diquark coupling (and on the formfactor of the mo¬ 
mentum space regularization, see 0 )- For all values of 
Gd no stable hybrid stars with a CFL phase have been 
found yet, see |T^|, and Refs, therein. We will restrict us 
here to the discussion of 2SC and NQ phases. 

The 2SC phase occurs at lower baryon densities than 
the CFL phase [H0- For applications to compact stars 
the omission of the strange quark flavor is justified by the 
fact that chemical potentials in central parts of the stars 
barely reach the threshold value at which the mass gap 
for strange quarks breaks down and they may appear in 
the system [18l |. 

It has been shown in 0 that a nonlocal chiral quark 
model with the Gaussian formfactor ansatz leads to an 
early onset of the deconfinement transition so that hybrid 
stars with large quark matter cores fll} can be discussed. 

In describing the hadronic part of the hybrid star, as in 
Jj], we adopt the Argonne Vl8+5v + UIX* model for the 
EoS jjSjj, which is based on recent data for the nucleon- 
nucleon interaction with the inclusion of a parameterized 
three-body force and relativistic boost corrections. 

Actually we continue to adopt an analytic parameter¬ 
ization of this model by Heiselberg and Hjorth-Jensen 
0, where the fit is done for n < 4 no with no = 
0.16 fm -3 being the nuclear saturation density. This 
EoS fits the symmetry energy to the original Argonne 
U18 + 5v + UIX* model in the mentioned density inter¬ 
val and smoothly incorporates causality constraints at 
high densities. The threshold density for the DU process 
is nj? u — 5.19 no, i.e. it occurs in stars with masses 
exceeding M,P U ~ 1.839 -M©). 


B. Cooling 

For the calculation of the cooling of the hadronic part 
of the hybrid star we use the same model as in The 
main processes are the medium modified Urea (MMU) 
and the pair breaking and formation (PDF) processes for 
our adopted EoS of hadronic matter. For a recent, more 
detailed discussion of these processes and the role of the 
3P 2 gap, see 0. 

The possibilities of pion condensation and of other so 
called exotic processes are suppressed since in the model 
|Z] these processes may occur only for neutron star masses 
exceeding M^ uark = 1.214 Mg. The DU process is irrel¬ 
evant in this model up to very large neutron star masses 
M > 1.839 Mq. The l<Sb neutron and proton gaps are 
taken the same as those shown by thick lines in Fig. 5 of 
Ref. @]. We pay particular attention to the fact that the 
3P 2 neutron gap is additionally suppressed by the factor 
0.1 compared to that shown in Fig. 5 of (zj. This sup¬ 
pression is motivated by the result of the recent work in 


>2J] and is required to fit the cooling data. 

For the calculation of the cooling of the quark core 
in the hybrid star we use the model introduced in |7]. 
We incorporate the most efficient processes: the quark 
direct Urea (QDU) processes on unpaired quarks, the 
quark modified Urea (QMU), the quark bremsstrahlung 
(QB), the electron bremsstrahlung (EB), and the mas¬ 
sive gluon-photon decay (see [US])- Following f26) we 
include the emissivity of the quark pair formation and 
breaking (QPFB) processes. The specific heat incorpo¬ 
rates the quark contribution, the electron contribution 
and the massles and massive gluon-photon contributions. 
The heat conductivity contains quark, electron and gluon 
terms. 

The 2SC phase has one unpaired color of quarks (say 
blue) for which the very effective quark DLI process works 
and leads to a too fast cooling of the hybrid star in dis¬ 
agreement with the data We have suggested to as¬ 
sume a weak pairing channel which could lead to a small 
residual pairing of the hitherto unpaired blue quarks. We 
call the resulting gap Ax and show that for a density de¬ 
pendent ansatz 


A x = A 0 exp 



(^) 


(1) 


with /r being the quark chemical potential, /i c = 330 
MeV. Here we use different values of a and Aq, which are 
given in Table 1 characterizing the model. 

The physical origin of the X-gap remains to be identi¬ 
fied. It could occur, e.g., due to quantum fluctuations of 
color neutral quark sextett complexes [22j]. Such calcu¬ 
lations have not yet been performed with the relativistic 
chiral quark models. 

For sufficiently small Gd, the 2SC pairing may be in¬ 
hibited at all. In this case, due to the absence of this 
competing spin-0 phase with large gaps, one may invoke 
a spin-1 pairing channel in order to avoid the DU prob¬ 
lem. In particular the color-spin-locking (CSL) phase jjgj 
may be in accordance with cooling phenomenology as all 
quark species are paired and the smallest gap channel 
may have a behavior similar to Ecp m , see [22[ • A consis¬ 
tent cooling calculation for this phase, however, requires 
the evaluation of neutrino emissivities (see, e.g. |2£| and 
references therein) and transport coefficients, which is 
still to be performed. 

Gapless superconducting phases can occur when the 
diquark coupling parameter is small so that the pairing 
gap is of the order of the asymmetry in the chemical po¬ 
tentials of the quark species to be paired. Interesting 
implications for the cooling of gapless CFL quark mat¬ 
ter have been conjectured due to the particular behavior 
of the specific heat and neutrino emissivities 0. For 
reasonable values of Gd, however, these phases do oc¬ 
cur only at too high temperatures to be relevant for late 
cooling, if a stable hybrid configuration with these phases 
could be achieved at all 0- 
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FIG. 1: Hybrid star cooling curves for Model I. Different lines 
correspond to compact star mass values indicated in the leg¬ 
end (in units of Mg), data points with error bars are taken 
from Ref. J|. For the explanation of shaded areas, see text. 



FIG. 2: Same as Fig. Qfor Model II. 


The weak pairing channels are characterized by gaps 
typically in the interval 10 keV -Fl MeV, see discussion 
of different attractive interaction channels in paper 
In the figures we present T-t plots for four models used 
in this paper. On each plot data points for known cooling 
NSs are added (see details in {§]). The hatched trapeze¬ 
like region represents the brightness constraint (BC). For 
each model nine cooling curves are shown for configu- 


FIG. 3: Same as Fig. □ Model III. 



FIG. 4: Same as Fig. |l]for Model IV. 

rations with mass values corresponding to the binning 
of the population synthesis calculations explained in the 
next section. 

Clearly, all models satisfy the BC. As for the T-t test 
the situation is different. The Model II does not pass the 
test because even the highest mass configuration (which 
corresponds to the coolest HyS) cannot explain the lowest 
data points. So, in the Table 1 it is marked that the 
model does not pass the test. 

In this work we want to introduce a more detailed mea¬ 
sure for the ability of a cooling model to describe obser- 
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vational data in the temperature-age diagram. We assign 
five grey values to regions of compact star masses in the 
T — t diagram which encode the likelihood that stars in 
that mass interval can be found in the solar neighbor¬ 
hood, according to the population synthesis scenario, see 
Fig. 0 The darkest grey value, for example, corresponds 
to the mass interval 1.35 -=-1.45 M 0 for which the popu¬ 
lation sysnthesis predicts the most objects. According to 
this refined mass spectrum criterion a cooling model is 
optimal when the darkness of the grey value is in accor¬ 
dance with the number of observed objects in that region 
of the temperature-age diagram. This criterion is ideally 
fulfilled for Model IV, where with only one exception all 
objects are found in the two bands with darker grey val¬ 
ues whereas for Models I - III about half of the objects 
are situated in light grey or even white regions. 


III. POPULATION SYNTHESIS SCENARIO 

Population synthesis is a frequently used technique in 
astrophysics described, e.g., in the review |33] where fur¬ 
ther references can be found. The idea is to construct an 
evolutionary scenario for an artificial population of cer¬ 
tain astronomical objects. The comparison with observa¬ 
tions gives the opportunity to test our understanding of 
evolutionary laws and initial conditions for these sources. 

The scenario that we use in this paper is nearly identi¬ 
cal to the one used in Paper I. We just briefly recall the 
main elements and then describe the only small difference 
in the mass spectrum. 

The main ingredients of the population synthesis 
model we use are: the initial distribution of NSs and 
their birth rate; the velocity distribution of NSs; the mass 
spectrum of NSs; cooling curves and interstellar absorp¬ 
tion. 

In this series of papers in which we use the population 
synthesis model as a test of the theory of thermal evo¬ 
lution of NSs, we assume that the set of cooling curves 
is the most undetermined igredient. So, we make an at¬ 
tempt to test it. The cooling curves used in this paper 
are described in Sec. 2. 

We assume that NSs are born in the Galactic disc and 
in the Gould Belt. The disc region is calculated up to 3 
kpc from the Sun, and is assumed to be of zero thickness. 
The birth rate in the disc part of the distribution is taken 
as 250 NS per Myr. The Gould Belt is modeled as a flat 
disc-like structure with a hole in the center (see the Belt 
description in [34j). The inclination of the Belt relative 
to the galactic plane is 18°. The NS birth rate in the 
Belt is 20 per Myr. 

The velocity distribution of NSs is not well known. In 
our calculations we use the one proposed by |3fj| . This 
is a bimodal distribution with two maxima at ~ 127 and 
~ 707 km s _1 . Recent results question this bimodality 
136). However, since the time scales in our calculations 
typically are not very long, the exact form of the distri¬ 
bution is not very important. 


For the calculation of the column density towards a 
given NS we use the same approximation as we used be¬ 
fore. It depends only on the distance to the galactic cen¬ 
ter and the height above the galactic plane (see Fig. 1 in 
H3). The detailed structure of the interstellar medium 
(ISM) is not taken into account, except the Local Bub¬ 
ble, which is modeled as a 140-pc sphere centered on the 
Sun. 

The mass spectrum is a crucial ingredient of the sce¬ 
nario. The main ideology of its derivation is the same as 
given in [38l |. At first we take all massive stars which can 
produce a NS (spectral classes B2-08) from the HIPPAR- 
COS catalogue with parallaxies < 0.002 arcsec. Then 
for each spectral class we assign a mass interval. In the 
next step using calculations by [3?)l | we obtain the baryon 
masses of compact remnants. Then we have to calculate 
the gravitational mass. In this paper, unlike our pre¬ 
vious studies where we just used the formula from jT(|| 
Afbar — Afg rav = a Mg rav with a = 0.075, we use M grav 
accurately calculated for the chosen configuration. Since 
the approximation from ,4Q] is very good the difference 
with the mass spectrum we used in Paper I is tiny, and 
appears only in three most massive bins of our spectrum. 

In this paper we slightly rebinned the mass spectrum in 
order to have a better coverage of the cooling behaviour 
for the chosen configurations. As before we use eight 
mass bins defined by their borders: 1.05; 1.13; 1.22; 1.28; 
1.35; 1.45; 1.55; 1.65; 1.75 M 0 , see Fig. 0 The crit¬ 
ical mass for the formation of a quark core is close to 
1.22 Mq. Therefore, bins are chosen such that the first 
two represent purely hadronic stars. Bins are of different 
width. Outermost bins have a width 0.1 M 0 . We do 
not expect HyS with masses M < 1.05 M 0 . The upper 
boundary of the eighth bin lies close to the maximum 
mass allowed by the chosen configuration. 

Following the suggestion by [40] we make runs for two 
modifications of the mass spectrum. Except the usage 
of the full range of masses we produce, in addition, cal¬ 
culations for the truncated spectrum. In this case con¬ 
tributions of the first two bins are added to the third 
bin. This situation reflects the possibility that stars with 
M ~ 11 Mq can produce NSs of similar masses close to 
~1.27 Mq. 

As in Paper I we neglect effects of a NS atmosphere, 
and use pure blackbody spectra. As we do not address 
particular sources this seems to be a valid approximation. 

The population synthesis code calculates spatial tra¬ 
jectories of NSs with the time step 10 4 yrs. For each 
point from the set of cooling curves we have the surface 
temperature of the NS. Calculations for an individual 
track are stopped at the age when the hottest NSs (for 
all five models here this is a star with M = 1.1 Mq, unless 
the truncated mass spectrum is used) reaches the tem¬ 
perature 10 5 K. Such a low temperature is beyond the 
registration limit for ROSAT even for a very short dis¬ 
tance from the observer. With the known distance from 
the Sun and the ISM distribution we calculate the col¬ 
umn density. Finally, count rates are calculated using the 
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Mass 

FIG. 5: The adopted mass spectrum, binned over eight in¬ 
tervals of different widths. The non-truncated spectrum is 
shown (see text). 

ROSAT response matrix. Results are summarized along 
each individual trajectory. We calculate 5,000 tracks for 
each model. Each track is applied to all eight cooling 
curves. With a typical cooling timescale of about 1 Myr, 
we obtain ~ 4 • 10 6 “sources”. The results are then nor¬ 
malized to the chosen NS formation rate (290 NSs in the 
whole region of the problem). 

IV. NUMERICAL RESULTS 

In this section we present results of our calculations for 
four models, characterized by different sets of the two- 
parameter ansatz for the X-gap, Ecp d see Table 1. 
Log N-Log S curves are given for two values of the Gould 
Belt radius (Rbeit = 300 and 500 pc), and for two variants 
of the mass spectrum (full and truncated). 

The modeled Log N-Log S curves are confronted with 
data for close-by, young cooling NSs observed by ROSAT. 
This data set includes the Magnificent Seven (seven dim 
radio-quiet NSs), radio pulsars, Geminga and a geminga- 
like source (see the list and details in |41j). The error bars 
correspond to poissonian errors (i.e. square root of the 
number of sources). An important upper limit is added 
)42j which represents an estimate of unidentified cooling 
NSs in the ROSAT Bright Source Catalogue (BSC). 

Model I is the best model from Paper II. An important 
feature of this model is that cooling curves cover data 
points in the T-t plot very uniformly. 

The parameters of the Model II were specifically cho¬ 
sen in such a way, that it is possible to demonstrate the 
fact that the Log N-Log S test can be successful for a set 
that fails to pass the T-t test. Even for the highest pos¬ 
sible mass it is imposible to explain cold stars, but as the 


Log N-Log S test is not sensitive to what is happening 
with massive NSs it does not influence the results of the 
population synthesis. 

Model III is an attempt of a compromise to fulfill at 
least marginally all three tests. It has s smaller gap than 
the Model I, and unlike Model II it has a non-zero value 
of a. It can explain all data points in the T-t plot within 
an available mass range resorting, however, to the very 
unlikely objects with masses above 1.5 M 0 . 

Models IV assumes a steeper density dependence of 
the X-gap than all previous ones. It is thus possible to 
spread the set of cooling curves over the existing cooling 
data already for mass variations within the range of most 
probable mass values 1.25 ± 0.25 M 0 in Fig. 0 Using 
these models it is possible to describe even the Vela pulsar 
being a young, nearby and rather cool object, within this 
mass range. 


V. DISCUSSION 

In the present paper we were to find a set of parame¬ 
ters for which all three tests (T-t, Log N-Log S, and BC) 
can be successfully passed? Are these tests sufficient to 
constrain a model or is it necessary to assume additional 
constraints? 

For the T-t test it is necessary to cover the observed 
points by curves from a relatively wide range of masses 
which are consistent with known data on inital NS mass 
distribution (for example, data on masses of not-accreted, 
i.e. mainly secondary, companions of double NS sys¬ 
tems). For the Model III (see Fig. II, for example, this is 
not the case as a number of data points seems to corre¬ 
spond to a narrow mass range slightly below the critical 
mass 1.22 M 0 . On the one hand, all data points on Fig. El 
can be covered by cooling curves from the standard mass 
range (~ 1 - 1.5, M@). On the other hand, the interme¬ 
diate region (logT ~ 6 - 6.2, logf ~ 3 - 4) corresponds 
to a narrow mass range: ~ 1.22 - 1.26 M 0 . This is 
not a dramatic disadvantage, especially if the hypoth¬ 
esis discussed in piiil. |4Q] that stars with masses below 
~ 11 Mq form remnants of nearly the same mass close 
to the range above. Still, this property of cooling curves 
should be mentioned. Model IV (Fig.QJ gives a more ap¬ 
propriate description from the point of view of the mass 
distribution. 

For the case of the Log N-Log S test we’d like to note 
that the value Rbeit=300 pc is more reliable. So, the 
Model I, for which at bright fluxes we see an overpredic¬ 
tion of sources for this value of the Gould Belt radius, 
can be considered only as marginally passing the test. 
Other models do better. Especially, models III and IV. 
For them Log N-Log S curves for Rb e it = 300 pc match 
well the data points, leaving a room for few new possible 
identifications of close-by cooling NSs (active search is 
going on by different groups in France [J3, m , in Ger¬ 
many 4T»|. in Italy 16| and in the US.\I li1| l. 

As in our Log N-Log S calculations we use a particu- 
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FIG. 6: LogN-LogS distribution for Model I. Four variants 
are shown: Rbeit = 500 pc and truncated mass spectrum 
(full line), Rbeit = 500 pc and non-truncated mass spectrum 
(dotted line), Rbeit = 300 pc and truncated mass spectrum 
(dash-dotted line), and finally Rbeit = 300 pc (dashed line) 
for non-truncated mass distribution. 



S [cts s" 1 ] 


FIG. 7: LogN-LogS distribution for Model II. Line styles as 
in Fig. El 


lar mass spectrum in which there are nearly no objects 
with M ~ 1.4 - 1.5 M 0 , we have a strong constraint 
on properties of a set of cooling curves which can satisfy 
all three tests. The position of the critical curve which 
devides hadronic stars from HyS is fixed by the chosen 



FIG. 8: LogN-LogS distribution for Model III. Line styles as 
in Fig. El 
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FIG. 9: LogN-LogS distribution for Model IV. Line styles as 
in Fig. El 


configuration. If curves for masses lip to 1.4 M e lie too 
close to the critical one, then we overpredict the num¬ 
ber of sources on the LogN-LogS plot. If oppositely, we 
move curves for 1.3 -1.4 M 0 down (this is achieved by 
increasing the parameter a and thus making the density 
dependence of the X-gap steeper) then a narrow range of 
masses becomes responsible for a wide region in the T -1 
diagram. 

A solution could be in changing the exponential de- 
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pendence in Eq. 1 to the power-law. We plan to study 
this possibility in future. 

The combined usage of all three tests can put addi¬ 
tional constraints on the mass spectrum of compact ob¬ 
jects. For example, if we look at Fig. 0 it is clear that 
small masses (M < 1.2 M 0 ) are necessary to explain hot 
objects with ages ~ 10 3 - 10 4 yrs. If in such a case only 
a model with a truncated mass spectrum is able to ex¬ 
plain the LogN-LogS distribution of close-by NSs then 
the model is in trouble. On other hand, if an explanation 
of the Log N-Log S deserves stars from low-mass bins, but 
cooling curves for these stars are in conradiction with the 
BC, then, again, the model has to be rejected. 

In this series of our studies of the local population 
of cooling compact objects we use the mass spectrum 
which should fit to the solar neighbourhood enriched 
with stars from the mass range 8-15 M 0 . Mass spec¬ 
trum of all galactic newborn NSs can be different, but 
not dramatically. The number of low-mass stars (~ 1 - 
1.3 Mq) can be slightly decreased in favour of more mas¬ 
sive stars. However, compact objects with M ~ 1.5 M 0 
anyway should significantly outnumber more massive ob¬ 
jects. This claim has some observational support. 

Unfortunately, a mass determination with high preci¬ 
sion is available only for NSs in binary systems. Com¬ 
pact objects in X-ray binaries could accrete a significant 
amount of matter. Also mass determinations for them 
are much less precise than for radio pulsar systems. So, 
we concentrate on the latter. For some of the radio pul¬ 
sars observed in binaries, accretion also played an impor¬ 
tant role. Without any doubts masses of millisecond pul¬ 
sars do not represent their initial values. However, there 
is a small number of NSs with well determined masses, 
for which it is highly possible that these masses did not 
change significantly since these NSs were born (data on 
NS masses can be found, for example, in HI [Mi] and 
references therein). These are secondary (younger) com¬ 
ponents of double NS systems. According to standard 
evolutionary scenarios these compact objects never ac¬ 
creted a significant amount of mass (as when they formed 
the primary component already was a NS). Their masses 
lie in the narrow range 1.18 - 1.39 Mq. Primary com¬ 
ponents of double NS binaries could accrete during their 
lifetime. However, this amount of accreted matter can¬ 
not be large as these are all high-mass binaries. Masses 
of these NSs are all below 1.45 Mq. This is also an im¬ 
portant argument in favour of the statement that initial 
masses of most NSs are below ~ 1.4 - 1.5 Mq. 

The recently discovered highly relativistic binary pul¬ 
sar .11906+0746 HI is the nineth example of such a sys¬ 
tem. The total mass is determined to be 2.61 ± 0.02 Mq. 
The pulsar itself is a young, not millisecond, object. It 
should not increase its mass due to accretion. So, we can 
assume that it is at least not heavier than the second - 
non-pulsar - component (we neglect here the possibility 
that the companion is a massive white dwarf, still this is 
a possibility), than we obtain that its mass is ~ 1.3 Mq. 
Nine examples (without any counter-examples) is a very 


good evidence in favour of the mass spectrum used in our 
calculations. Of course, some effects of binary evolution 
can be important (for example, the mass of the core of a 
massive star can be influenced if the star looses part of its 
mass due to mass transfer to the primary companion or 
due to common envelope formation), and so for isolated 
stars (or stars in very wide binaries) the situation can 
be slightly different. However, with these observational 
estimates of initial masses of NSs we feel more confident 
using the spectrum with a small number of NSs with 
M ~ 1.4 - 1.5 Mq. 

Brighter sources are easier to discover. So, among 
known cooling NSs the fraction of NSs with masses 
1 Mq ~ M ~ 1.5 Mq should be even higher than in the 
original mass spectrum. So, we have the impression that 
it is necessary to try to explain even cold (may be with an 
exception of 1-2 coldest) sources with M ~ 1.4 - 1.5 Mq. 
Especially, the Magnificent seven and other young close- 
by compact objects should be explained as most typical 
representatives of the whole NS population H3|- We want 
to underline that, even being selected by their observ¬ 
ability in soft X-rays, these sources form one of the most 
uniform samples of young isolated NSs. In this sense, 
the situation as in Fig. Q] where a significant number of 
sources are explained by cooling curves corresponding to 
1.5 Mq ~ M ~ 1.7 Mq should be considered as a disad¬ 
vantage of the model. Particularly, Vela, Geminga and 
RX .11856-3754 should not be explained as massive NSs. 
These all are young close-by sources, and the probability 
that so near-by we observe young NSs which come out of 
few percent of the most massive objects of this class is 
low. 

All the above gives us the opportunity to formulate 
the conjecture of a mass spectrum constraint: data points 
should be explained mostly by NSs with typical masses. 
For all known data this typical means 1.1 M 0 ~ M ~ 
1.5 M 0 . 

The LogN-LogS test is the only one that takes into 
account the mass spectrum explicitely. This is an addi¬ 
tional argument in favour of using this test together with 
others. 

Taking all together, we conclude that Model IV is the 
best among studied examples from the point of view of 
all three tests plus mass constraint. 


VI. CONCLUSION 

We made a preliminary study of cooling curves for HyS 
based on the approach in Paper II which suggests that if 
quark matter occurs in a compact star, it has to be in the 
2SC+X phase, where the hypothetical X-gap still lacks a 
microscopical explanation. All three tests of the cooling 
behavior (LogN-LogS, T-t, BC) are applied. Four mod¬ 
els defined by the two-parameter ansatz for the X-gap 
were calculated. Model II with a density-independent X- 
gap could directly be excluded since it was not able to 


explain some cooling data including Vela at all. Two of 
the models (I and III) successfully passed two tests and 
marginally the third. None of these models could ex¬ 
plain explain the temperature-age value for Vela within 
a typical mass range, i.e. for a Vela mass below 1.5 M®, 
However, with a steeper density dependence of the X- 
gap than suggested in Paper II, we were able to fulfill 
all 4 constraints and exemplified this for model IV. To 
conclude, HyS with a 2SC+X quark matter core appear 
to be good candidates to explain the cooling behaviour 
of compact objects, although a consistent theoretical ex¬ 
planation of the hypothetical X-gap and its steep density 
dependence has still to be developed. 
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Model 

Ao [MeV] 

a 

BC 

T - t 

Log N - Log S 

Mtyp < 1.5 M© 

All tests 

I 

1 

10 

+ 

+ 

o 

- 

- 

II 

0.1 

0 

+ 

- 

+ 

- 

- 

III 

0.1 

2 

+ 

o 

+ 

- 

- 

IV 

5 

25 

+ 

+ 

+ 

+ 

+ 


TABLE I: Parameter values Ao and a for Models I - IV defining the density dependence of the X-gap introduced in Eq. 0. The 
symbols +, - and o denote whether or not a given model passes the proposed cooling tests, fails or is indifferent, respectively. 














